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Abstract 

Closed form, analytical results for the finite-temperature one-body density matrix, and Wigner function 
of a d-dimensional, harmonically trapped gas of particles obeying exclusion statistics are presented. As 
an application of our general expressions, we consider the intermediate particle statistics arising from the 
Gentile statistics, and compare its thermodynamic properties to the Haldane fractional exclusion statistics. 
At low temperatures, the thermodynamic quantities derived from both distributions are shown to be in 
excellent agreement. As the temperature is increased, the Gentile distribution continues to provide a good 
description of the system, with deviations only arising well outside of the degenerate regime. Our results 
illustrate that the exceedingly simple functional form of the Gentile distribution is an excellent alternative 
to the generally only implicit form of the Haldane distribution at low temperatures. 
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I. INTRODUCTION 

In his now famous 1991 paper, F. D. M. Haldane | 
Pauli exclusion principle, that leads to particle statistics which continuously interpolate between 
the Bose and Fermi statistics. In the Haldane fractional exclusion statistics (FES), one constructs 
a generalized exclusion principle through the single particle dimension of the iV-th particle, in the 
presence of N — 1 other identical particles, viz., M\ 

d 9 N = d-g(N-l), (1) 

where we note that g = corresponds to bosons, and g = 1 to fermions. The constant g is the 
so-called FES of a particle, and by definition is given by 

In Eq. (2), Ad denotes the change in the dimension of the single particle space, and AiV is the 
change in the number of particles, with the proviso that the size and boundary conditions of the 
system are unchanged. The FES parameter, g, is then a measure of partial Pauli blocking, and can 
quite generally take on arbitrary values g > 0, although we will focus on < g < 1. Furthermore, 
FES is a consequence of state counting arguments made in Hilbert space, and so is valid for 
arbitrary spatial dimensions, yfl 

Following Haldane's work, R. Ramanathan [4], Dasineres de Veigy and Ouvry [5], Wu [6] and 
Isakov tn\ independently examined the thermodynamic properties of an ideal FES gas, and derived 
the now well-known result for the average occupancy of a gas of particles obeying ideal FES in 
the grand canonical ensemble, viz., 

where T is the temperature, q is the single particle energy, k B is the Boltzmann constant, and \i is 
the chemical potential. The function /(e (e4_At )/ fcsT ) is determined by (x = e^ ei ~^^ kBT ) 

f(xy(i + f(x)y-v = x. (4) 

It is easily seen that for g — 0, we recover the usual ideal Bose distribution, while for g = 1, we 
obtain the Fermi distribution. 

In general, closed form expressions for n(ej) in the Haldane FES, Eq. (3), are not possible; that 
is, n(ej) is generally only given implicitly through Equation (4). However, for the special case 



of g = p/q where p and q are co-prime, explicit expressions may be obtained. [8] Regrettably, 
such distributions quickly become difficult to work with both numerically and analytically. For 
example, with g = 1/2 (semions), and g — 1/3, one obtains Q8D 

2 

nUi) = . , (5) 

and 

n(e * } = fe+A-i ' (6) 

h=[2y/y(y-l) + 2y-l] 1 / 3 , y = 2 +1, (7) 
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respectively. It is therefore reasonable to investigate if there is an alternative distribution to Eq. (3), 
which is able to capture the essential thermodynamic properties of the Haldane FES, while still 
possessing the desirable property of having an explicit functional form. 



Remarkably, over 70 years ago, the pioneering work of G. Gentile |^] on intermediate parti- 
cle statistics may have already provided a possible answer to this question. Gentile's work was 
founded on a quantum phase-space approach, in which a single quantum cell may accommodate 
up to 1/ a w particles with the same energy. Utilizing the method of most probable distribution, 
Gentile obtained what we shall refer to as the Gentile exclusion statistics (GES) distribution, viz., 

nMff.) = : ™ (X) 

where w = 1 ~^* w . It is clear that for a w = and a w — 1, Eq. (8) exactly reproduces the standard 
Bose and Fermi distributions, respectively. More importantly, Eq. (8) also provides a simple, form 
invariant, expression for the average occupation number, n^ aw \ei), V < a w < 1. We wish to 
point out that Eq. (8) appears to have been recently "rediscovered" by Q. A. Wang et. al Jlpj. who 
based their analysis on the grand partition function, also supposing n max = l/a w for the maximum 
occupation number. In fact, the distribution derived in Ref. Jioll is identical to what is obtained in 
Gentile statistics, Eq. (8), although Wang et al. do not appear to be aware of this fact. 

One of the objectives of this paper is to investigate the viability of the GES explicit distribution, 
viz., Eq. (8), as an alternative to the implicit FES distribution given by Equations (3) and (4). We 
note that this is a meaningful comparison, since both GES and FES are rooted in the generaliza- 



tion of the Pauli principle, with a well defined method for the counting of states. [|13|1 In order to 
facilitate this goal, we will focus our attention to a (/-dimensional gas of ideal particles obeying 
arbitrary statistics at finite-temperature, confined to a harmonic oscillator trap, V(R) = mu 2 R 2 /2, 



with R = \Jxf + x\ + • ■ ■ + x\ being the d-dimensional hyper-radius. The motivation for study- 
ing this system lies in its possible connection to current experiments on harmonically trapped, 
ultra-cold quantum gases, along with the models relatively simple analytical properties. 

To this end, the rest of our paper is organized as follows. In Sec. II, we will present finite 
temperature, closed form analytical expressions for the d-dimensional one-body density matrix 
(ODM) and Wigner function obeying general exclusion statistics. Then, in Sec. Ill, we make use 
of the Wigner function to construct a variety of thermodynamic properties without restriction to 
any specific statistics. In Sec. IV, we narrow our focus to FES and GES, so that we may make a 
detailed comparison of these two distributions at finite, and zero-temperature. In Sec. V we present 
our concluding remarks. 



II. FINITE TEMPERATURE ONE-BODY DENSITY MATRIX AND WIGNER FUNCTION 

In this section, we will provide closed form expressions for the finite-temperature ODM and 
Wigner function obeying exclusion statistics, in arbitrary dimensions for a harmonically trapped 
gas. In what follows, we will denote the general exclusion statistics parameter by a, such that 
< a < 1 defines the fermionic sector. This is in fact a generic feature of any exclusion statistics 
distribution, which is required to continuously interpolate between Bose and Fermi statistics (see 
also Eq. (15) below). The finite-temperature Wigner function is subsequently used to evaluate 
a variety of thermodynamic properties, such as the spatial density, momentum density, kinetic 
energy density, and form factor. 



A. One-body density matrix 



The d-dimensional, finite-temperature ODM, s;T), for a system obeying arbitrary 

statistics is obtained by taking the two-sided inverse Laplace transform (ILT) of the finite- 
temperature Bloch-density matrix, Ill4ll 



.(4 



(9) 



where 



4 d) (r, S ;/3) = CT(r,s;/3)x 



H^(/3;T) 



(10) 
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and j3 is the generally complex variable conjugate to fi. We have also introduced the center-of-mass 
and relative coordinates 

and s = ri — r 2 , (11) 
h — m — co — 1. The factor, 



respectively. Hereby, we shall use units such that k B = ti = m = u 
H^ a \(3;T)/ (3, in Eq. (10) is a thermal weighting factor, taking into account the statistics, and 
C^ d) (r, s; (3) is the normal (T = 0) Bloch-density matrix. |l4] 

The ILT in Eq. (9) will be a convolution between the /3-dependence of Cq (r, s; f3), and 
H( a \(3;T)/(3. Therefore, it would appear that the thermal factor must be known explicitly in 
order to obtain the ODM of the system. However, we now point out the following defining prop- 
erty of the two-sided ILT for any given thermal factor (f3;T)/ (3, which we write as 



ST, 1 



(12) 



InEq. (12), n^(e) is the distribution associated with the statistics of the particles. For example, 
let us consider the Bose and Fermi distributions 



1 



and 



n B [e) 



n F le 



1 



e (e-^)/T + J 



respectively. Generalizing the above cases, one may write Ill5h 

1 



n {a) (e) 



(13) 



(14) 



(15) 



e (e-R)/T + a ' 

so that Bose (a = —1), Boltzmann (a = 0), and Fermi (a = 1) statisics are all represented by a 
universal distribution. Focusing on the Bose and Fermi distributions, explicit expressions for the 
thermal factors H^ a \(3; T)/ (3 are known; namely, for bosons jlq] 

H B ((3;T) ttT 



whereas for fermions, 



(3 tan(vr/3T) 
H F ((3;T) ttT 







ain(7r/3T) 



It can be readily shown by direct calculation that Ill6h 



B -i 



(3 



1 



exp(— /i/T) — 1 



n B (0) 



(16) 



(17) 



(18) 



and 



B7, 1 



H F {I3-T) 



1 



n F (0) 



(19) 



exp(—fi/T) + 1 

Therefore, once the /3-dependence of Cq (r, s;/3) is known, an application of the convolution 
theorem for Laplace transforms will, at least in principle, be able to provide us with the finite- 
temperature ODM via Equation (9). Note that the convolution integral may still be very difficult 
to evaluate analytically if the /3-dependence coming from the normal Bloch-density matrix is com- 
plicated. Indeed, the analytical evaluation of the ILT may only be feasible for specific dimensions. 

However, we now make a critical observation. If Cg (r, s; (3) can be written such that the /in- 
dependence is exponential, the shifting property of the Laplace transform 111711 may be used to find 
a universal expression for the finite-temperature ODM, which is unchanged by the dimension or 
statistics under consideration. To wit, we note that by the shift property, we have 



n (a) (7) 



(20) 



where 7 is real and positive. In other words, if the Bloch density matrix is purely exponential in its 
/3-dependence, the ODM may be found for any statistics without requiring an explicit expression 
for the thermal factor, H^ a \(3; T) / (3. It is then highly desirable to try to express the /3-dependence 
of Cq^ (r, s; /3) as a pure exponential. This goal is actually achievable for the case of rf-dimensional 
harmonic confinement, where we obtain the following expression for the normal Bloch-density 
matrix, viz., UM 



cf (r, S ;/3) 



9s 

n d/2 



OO CXD 



EE(- 1 ) nL n /2 " 1 ( 2 ^) L f" 1 ( s2 /2)e- 



^+s^A)-{e n +k)fj 



(21) 



n=0 k=0 

where e n = (n + d/2) is the <i-dimensional spectrum of an isotropic harmonic oscillator potential, 
and plays the role of 7 in Equation (20). The quantity, g s in Eq. (21) denotes the spin degeneracy, 
and L a n (x) are the associated Laguerre polynomials. QJ/ZU 

It then immediately follows that the ILT in Eq. (9) may be performed without requiring an 
explicit expression for the thermal factor, by using the general result 



B' 1 



n {a) (e r 



k). 



(22) 



In order to clarify, and illustrate the above analysis, let us again consider the Bose and Fermi 
statistics, from which Eqs. (18) and (19) provide us with the appropriate ILTs for the thermal 
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factors, if T) / /?. We may evaluate the ILT piece in Eq. (9) for bosons by brute force using 
the convolution theorem for ILTs, viz., 

1 



zr 1 

A* 



-(e„+fc)/3_ 



-7TT 



tan(7r/3T) 



rfr tf(r - (e B + fc)) 



exp[(r-/x)/T] - 1 



n B {e n + k) , 



(23) 



and similarly for fermions 



B -i 



-(e n +k)0_ 



TXT 



sm(irf3T) 



dr 5{t — (e n + k)) 



exp[(r - fjL)/T\ + 1 



e ( £n+fc _ M )/T _|_ j 

n F (e„ + fc) . 



(24) 



Notice that a direct application of Eq. (20) immediately leads to the same result, without explicit 
knowledge of H^(^; T)//3. 

We may then write down the general expression for the finite-temperature ODM of a harmoni- 
cally trapped gas, appropriate for general exclusion statistics, as 



oo oo 



P id \r,s;T) 



7T 



J2 ^(-l) n ^ /2 - 1 (2r 2 )Lf ~\s 2 /2)e-^ +s2 ^n^(e n + k) . (25) 

n=0 fc=0 



d/2 



Thus, for bosons, Eq. (25) would read exactly as above, but with n^(e n + k) — > n B (e n + k). 
Similarly, for fermions, Eq. (25) still holds, but with rS a \e n + k) — > np{^ n + k). We must 
emphasize that this universal expression for the finite-temperature ODM is only possible owing to 
the special decomposition of Cq (r, s; (3), such that the f3 dependence is strictly exponential. This 
is by no means a trivial result, and would be difficult, if not impossible to establish by starting 
with the single-particle harmonic oscillator eigenstates in c/-dimensions. In fact, for any other 
form of Cq (r, s; (3), the temperature dependence in Eq. (25) changes with dimensionality, and 
the universal representation of the ODM is lost, as illustrated in Reference II 1911 . In addition, note 
that in Eq. (25), both the center-of-mass, r, and relative coordinate, s, are treated on equal footing, 
resulting in a clean separation of the variables. This form for the ODM is useful for analytical 
calculations where separate integrations over r and s may need to be performed. 
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B. Wigner function 



The d-dimensional Wigner function, W^(r,p; T), may now be obtained via a Fourier trans- 
form of the ODM, Eq. (25), with respect to the relative coordinate. The Wigner function is a useful 
tool for the phase-space formulation of quantum mechanics, [20Q and as we shall see below, also 
simplifies analytical calculations for various thermodynamic properties of the system. Specifically, 
by definition, 

W (d) (r,p;T) = J d d s p {d \r,s;T)e- tp - s , (26) 

where d d x = 27r d ^ 2 x d ~ 1 /T[d/2). Given that the ODM only depends on the magnitude of the 
coordinates, all angular integrals may be immediately performed, thereby allowing us to write 
Eq. (26) as H 

POO / 1 \ ^/ 2— 1 

W {d \r,p;T) = (27r) d / 2 / ds p {d) (r, s;T)( — ) Jd,(sp)s d ~ 1 , (27) 

Jo \PsJ 2 



where J n (x) is a Bessel function of the first kind. Ill7h The integral in Eq. (27) has already been 



addressed in an earlier work, Q2JJ] and following the same analysis, we readily obtain the desired 
result 

oo oo 

W^(r,p;T) = 2^^^(-l)^^/ 2 - 1 (2r 2 )Lf" 1 (2p 2 )e-( r2+ ^)n^(e n + A ; ). (28) 

n=0 fc=0 

Similar to Eq. (25), there is once again a clean separation of the variables in Eq. (28), which in this 
case are the spatial and momentum variables. The utility of this form for the Wigner function will 
be illustrated below. 

The finite temperature expressions given by Eqs. (25) and (28) are valid for any dimensionality, 
any flavour of exclusion statistics, and represent the main analytical results of this paper. 

HI. FINITE TEMPERATURE RESULTS 

Here, we make use of the Wigner function developed above to construct several thermody- 
namic quantities of interest. The results presented here serve to generalize the Bose and Fermi 



expressions presented elsewhere in the literature. Q2JJ] 
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A. Spatial density 



The spatial density is obtained from the Wigner function via 

P {d \r-T) = JL J d d pW^(r,p;T) . (29) 



Using the integral Ill7h 



poo i 

/ = / dx Ll /2 -\2x 2 )e- x2 x d ~ 1 = 
Jo 2 

Eq. (29) evaluates to 



.T(k + d/2) 



(30) 



(k + d/2-l\ ., , 

n=0 fc=0 V ^ / 



pW(r;T) = ^^^(-ir [ ! )L^(2rV^)( en + fc). (31) 



Observe that the convenient separation of the r and p coordinates in Eq. (28) has allowed for an 
easy calculation of the spatial density. Of course, Eq. (31) may also be obtained by setting s = 
in Equation (25). 

B. Momentum density 

The finite temperature momentum density, IL^(p; T), is obtained by integrating over the coor- 
dinate variable, viz., 

n (d) (p;T) = JL J d d rW^(r,p;T) . (32) 
Once again, making use of Eq. (30), we obtain 



nM (p;T) = J^EB- 1 )* ( H + d/ n 2 1 ] Lf-\2p*)e-^(e n + k) . (33) 



TT d / 2 

n=0 k=0 



C. Kinetic energy density 

The finite temperature kinetic energy density, r^(r; T), is calculated according to 



r«(r; 



1 f°° ri 2 

T) = j2^rJ ddp wW ( r >n T) y • (34) 

Inserting the finite-temperature Wigner function, Eq. (28), into Eq. (34), and performing the inte- 
gration, leads to 
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tM( '' ;T > = ^7m^^ ( - 1) '' (d + 4 ' ) ^ /2 + t ' > L - /2 ''' 2r2 » e ' r, "'° ) '^ + ^ ' (35) 

V ' ' n=0 k=0 

From here, the kinetic energy can be obtained by integrating Eq. (1351) over all space, 

POO 

E km (T)= / d d rr^(r;T). (36) 
Jo 

Using tabulated integrals, |l7ll it is straightforward to show that 



= * ^ i^^(y^^ +t) . (37) 

^ / > n =0 k=0 



D. Form factor 



The form factor is simply the Fourier transform of the spatial density, and is given by, 

( g; T) — J d d r p {d) (r; T)e~ iq - r . (38) 
Following an identical analysis as for the evaluation of (r, p; T) yields Q2JJ] 



oo oo 



/» (g; T) = * £ £ ( ' + ^ ' ] L d { 2 ~ 1 (g 2 /2)e- 2 / 4 n(e n + /) 



n=0 (=0 



= ^t^gV^e^n^y. (39) 

m=0 

E. Zero temperature 

In the case of the fermionic branch, the T — > limit is easily obtained by noting that 

n^(e n + k)^ -Q(ef - (e n + fc)), (40) 

where the Fermi energy is now given by = yfa{M + d/2). Therefore, for any flavour of 
exclusion statistics, the zero-temperature ODM becomes 



n=0 fe=0 

^max femax 



a 



6( e W _ ( e „ + fc)) 



n=0 fc=0 



n=0 



^(-l)"^/ 2 - 1 (2r 2 )L£( S 2 /2)e-^ +s2 / 4 ) , (41) 
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where n max = \yfa(M + d/2)\ — d/2 and k„ 
T = Wigner function reduces to 



[y/a(M + d/2)\ - n - d/2. Similarly, the 



oo oo 



2\ 

a 



WW(r,p) = 2 d ( 7s ^^(-l)" +fc L^ 2 - 1 (2r 2 )Lf " 1 (2p 2 )e-( r - 2+ P 2 ) 

n=0 fc=0 

^max ^max 



e(4 a) 



^ + fc)) 



(42) 



n=0 fc=0 

All of the zero temperature results in the fermionic sector may be obtained from Equations (41) 
and (42). 



IV. APPLICATION 



While it may seem a little pedantic, we feel that it is useful to first illustrate the above discussion 
with a specific example, namely, the Gentile distribution, Equation (8). Even though we do not 
require an explicit expression for the thermal weighting factor, for the GES, it can readily be shown 
that 

HM(B:T) ttT 1 umT 

(43) 



(3 tan(vr/3T) u tan (uir/3T) ' 

With the above form for the GES thermal factor, one may then work out all of the two-sided ILTs 
explicitly, and readily confirm that this is equivalent to a direct application of Eq. (20), viz., 



h 1 



(3 



1 



1+Oi 



exp (*4£=fi) - 1 exp ((l±^) (^+j-M) ) _ i 



k). 



(44) 



Note that at T = 0, the fermionic branch {i.e., < a w < 1) of the GES becomes 



n 



f 



k)) 



(45) 



as mentioned above. We therefore have closed form, analytical expressions for the one-body 
density matrix, and Wigner Function for GES, given by Eq. (25) and (28), respectively, provided 
we take — > rS aw \ 

Let us continue the application of our results by also presenting the finite-temperature spatial 
density profiles (Eq. (31)) for a harmonically trapped system obeying GES, and comparing them 
to those obtained from the FES distribution. In this comparison, we identify a w with the FES 
parameter g. Our motivation is two-fold. First we wish to illustrate the quality of the much 
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simpler GES distribution when evaluated for a local quantity (as opposed to integrated thermody 



230). In 



namic quantities such as the chemical potential, specific heat or energy per particle 11221 . 
addition we would also like to examine how differences in the spatial densities may be used to 
probe the type of statistics exhibited by the system experimentally. In particular, we have in mind 
applications to ultra-cold, harmonically trapped Fermi systems in the unitary regime where there 
is suggestive evidence that the strongly interacting gas may be mapped to a noninteracting system 
obeying ideal FES. |24-E] 

In Fig. 1, we present the 3D spatial density profiles obtained from the GES (dashed curves) 
and FES (solid curves) at various temperatures, with g = a w — 1/2 (semions). Note that for 
this particular value of the statistical parameter, an explicit form for the distribution function is 
available, and is given by Equation (5). We note that at high temperature (panel (a)), the FES and 
GES spatial densities are in good agreement, with the only significant deviation occurring near 
the center of the trap. As the temperature is lowered, the agreement between the two densities 
improves. In the zero temperature limit, the two spatial densities are analytically identical, and 
the quantum mechanical shell oscillations become more prominent. As a reference, we have also 
included in the T = plot (panel (d)) the spatial density for the Fermi statistics (dotted line, a w = 
g = 1). We observe that the smaller statistical parameter also effectively serves to "bosonize" the 
particles, resulting in an increase in the density at the origin, a squeezing of the distribution in the 
tail region, and the diminished shell oscillations. 

In Fig. 2, we again present the spatial densities, but now with g = a w = 0.29. Our motivation 
for choosing this particular value of the statistical parameter lies in the earlier work of Bhaduri et 
al [|26|l in the context of a harmonically trapped ultra-cold Fermi gas, where it was argued that in 
the unitary regime, the strongly interacting system may be mapped onto a gas of particles obey- 
ing ideal FES. In their investigation, the statistical parameter was determined by fitting the finite 
temperature, theoretical FES energy per particle, E/N, and chemical potential to the experimen- 
tal data, thereby obtaining a "best fit" value of g = 0.29. It is important to note, however, that 
by fitting to only global quantities, local information is not included, which may be important in 
comparing to experimental data. We suggest that by also examining the local spatial density, one 
may be able to provide further evidence in support of the conjecture that the unitary Fermi gas 
obeys fractional statistics. 

Figure 2 once again illustrates that the GES and FES distributions are in good agreement, 
and as the temperature is lowered, the agreement improves. As in Fig. 1, the zero temperature 
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Figure 1: Three-dimensional, finite temperature spatial density profiles for N = 420 particles with g = 
a w = l/2 and g s = 2. (a) T = 3, (b) T = 2, (c) T = 1, (d) T = 0. Solid curves are FES, dashed curves are 
GES, and the dotted line in panel (d) is the T = Fermi density. Units are such that ks = h = m = cj = 1 
as discussed in the text. 



limit results in identical spatial densities between the two distributions. It is also clear that the 
reduced statistical parameter leads to more boson-like behaviour of the particles, as evidenced by 
the significant squeezing of the cloud, and the increased density in the central region of the trap, 
especially when contrasted with the Fermi density (panel (d), dotted line). It would therefore be 
interesting to examine the experimental density distribution of the trapped Fermi gas in the unitary 
regime, and compare it to the theoretical predictions of the GES and FES with g = a w = 0.29 
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Figure 2: As in Fig. 1. g = a w = 0.29. 



presented here. Owing to the similarities in the spatial distributions at low temperatures however 
(the regime of ultra-cold gases), it is unlikely that one would be able to determine the specific 
kind of the statistics obeyed by the particles. Nevertheless, the spatial density may well serve 
as another "smoking gun" signature that the system in the unitary regime is indeed exhibiting 
fractional statistics. 

We would like to further mention that the agreement between the FES and GES extends to all 
values of the statistical parameter, and is also evident in other thermodynamic quantities, which 
we have not detailed here. 
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V. CLOSING REMARKS AND CONCLUSIONS 



We have presented closed form, universal analytical expressions for the finite temperature one- 
body density matrix, and Wigner function of a rf-dimensional, harmonically trapped gas, obeying 
general exclusion statistics. These expressions, Eqs. (25) and (28), completely generalize results 
presented elsewhere, which were limited to Bose and Fermi statistics. 11191 . 12111 The universal forms 



of the one-body density matrix and Wigner function are only possible provided the normal Bloch- 
density matrix has its /3-dependence written in a purely exponential form. As a result, we have 
established that explicit knowledge of the thermal weighting factors, H^ a \/3;T)//3, previously 
thought to be necessary, [16, Jjl 21] for the evaluation of the one-body density matrix, are in fact 



not required. 

As an application of our results, we have examined the GES distribution, Eq. (8), which has re- 
cently been rediscovered, Jioll and proposed as an alternative to the more complicated distribution 
found independently by Ramanathan and others. (4h7|] Through an examination of the local spatial 
density at finite temperature, we were able to demonstrate that the GES is a good description of 
the harmonically confined gas obeying FES. Indeed, we have established that the low tempera- 
ture (ksT jfojj < 1) global and local thermodynamic properties derived from FES and GES are 
essentially indistinguishable. As a result, we note that it would be unlikely to experimentally as- 
certain the specific underlying fractional statistics of an ultra-cold Fermi gas in the unitary regime, 
as all such distributions will tend to lead to the same low temperature properties. In particular, 
any suggestion that the unitary Fermi gas obeys ideal FES is somewhat arbitrary, as almost iden- 
tical results will be found at low temperatures using some other exclusion statistics distribution 
which smoothly interpolates between Bose and Fermi statistics. For example, while we have not 
presented the details here, we have confirmed that the spatial densities obtained from the distri- 
bution given by Eq. (15) are indistinguishable from GES and FES at low temperatures, although 
noticeable differences from GES and FES do occur at higher temperatures. 

Given the excellent agreement between the GES and FES distributions at low temperatures, we 
conclude that the much simpler GES may be used with confidence in other studies where simple, 
analytical results for ultra-cold gases obeying FES are desired. 
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